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Abstract 

The  stability  of  a secondary  Tollmien-Schl ichting  wave,  whose 
wavenumber  and  frequency  are  nearly  one  half  of  those  of  a fundamental 
Tollmien-Schl ichting  instability  wave  (T-S  wave),  is  analyzed  by  using 
the  method  of  multiple  scales.  Under  these  conditions,  the  fundamental 
wave  acts  as  a parametric  excitor  for  the  secondary  wave.  When  the 
amplitude  of  the  fundamental  wave  is  small,  the  amplitude  of  the  secondary 
wave  deviates  slowly  from  its  unexcited  state.  However,  as  the  amplitude 
of  the  fundamental  wave  increases,  so  does  the  amplitude  of  the  secondary 
wave  even  in  the  regions  where  it  is  damped  in  the  unexcited  state. 
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I.  Introduction 

Secondary  instabilities  in  two-dimensional  boundary  layers  over 

flat  plates  is  analyzed.  One  of  the  many  disturbances  in  the  flow,  a 

Tollmien-Schlichting  instability  wave  (T-S  wave),  develops  and  starts  to 

amplify  at  a given  location  on  the  plate.  When  this  wave  grows  appreciably 

downstream,  the  flow  will  deviate  from  the  Blasius  flow  and  hence  it  may 

be  unstable.  When  this  occurs,  experimental  data  shows  that  the  higher 

harmonics  and  subharmonics  of  the  principal  T-S  wave  start  to  amplify 

faster  than  the  T-S  wave  itself.  It  is  therefore  predicted  that  tU 

higher  harmonics  and  the  subharmonics  of  the  T-S  wave  play  a big  role  in 

transition  to  turbulent  flows.  This  behavior  of  the  disturbances  is 

demonstrated  by  Kachanov,  Kozlov  and  Levchenko ^ and  Veasov,  Ginevsky  and 
2 

Karavosov  in  their  experimental  works  where  they  excited  a known 
distrubance  on  a flat  plate  by  a vibrating  ribbon  and  measured  the 
amplitudes  of  the  distrubance  spectrum  at  different  downstream  locations. 

Kachanov,  Kozlov  and  Levchenko^,  in  their  study  of  the  nonlinear 
development  of  a wave  in  boundary  layers,  divided  the  process  of  the 
development  of  a disturbance  Into  four  regimes:  (1)  a linear  and  slow 
nonlinear  development  of  dlstrubances  and  generation  of  higher  harmonics, 

(2)  an  Increase  In  the  amplitude  of  the  principal  wave,  a decrease  in 
the  amplitudes  of  its  harmonics  and  an  appearance  of  low  frequencies  and 
subharmonics  In  the  specturm,  (3)  an  Interaction  of  low-frequency  oscillations 
with  the  principal  wave  and  Its  higher  harmonics,  three-dimensional  and 
rapid  amplification  of  all  spectral  modes  and  (4)  a damping  of  singled 
out  harmonic  modes,  smoothing  of  the  spectrum  and  a transition  to  the 
turbulent  regime. 
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Klebanoff,  Tidstrom  and  Sargent  , in  their  experimental  study  of  the 
three  dimensional  nature  of  boundary- layer  instability,  observed  the 
secondary  instabilities  induced  by  nonlinear  interaction  between  a two- 
dimensional  fundamental  and  a three-dimensional  secondary  wave  at  a 
position  not  very  far  away  from  the  linear  region. 

In  the  present  study,  we  are  concerned  with  the  linear  and  two- 
dimensional  interaction  between  the  fundamental  T-S  wave  and  its  ^ 
subharmonic.  The  wavenumber  and  the  frequency  of  the  subharmonic  wave 
are  nearly  one  half  of  those  of  the  fundamental  wave.  The  initial  amp- 
litude of  the  fundamental  wave  is  varied  and  the  corresponding  amp- 
litudes of  the  subharmonic  wave  are  calculated  at  various  downstream 
locations.  It  is  found  that  the  amplitude  of  the  subharmonic  is  oscil- 
latory in  nature  and  is  not  amplified  appreciably  when  the  fundamental 
wave  is  considerably  small.  As  the  amplitude  of  the  T-S  wave  grows 

downstream,  its  subharmonic  grows  In  agreement  with  the  experimental 

4 

observations  of  Refs.  1 and  2.  von  Kerczek  analyzed  the  stability  of 
unsteady  boundary  layers  which  do  not  depend  on  the  streamwise  coordinate. 
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II.  Problem  Formulation 


l: 


We  consider  secondary  instabilities  of  a two-dimensional,  steady, 
incompressible  flow  past  a flat  plate.  The  equations  describing  the 
motion  are  the  unsteady,  dimensionless  Navier-Stokes  equations: 
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U -*■  U , V -*•  0 and  P -*■  p U2  as  y 
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a2  a2 

where  V2  = . Here,  x and  y are  made  dimensionless  by  using  a 


reference  length  6^,  the  time  is  made  dimensionless  by  using  and 


the  velocities  are  made  dimensionless  by  using  the  freestream  velocity 


U^.  The  Reynolds  number  is  defined  as  R = UjSp/v,  where  v is  the  fluid 


kinematic  viscosity. 

A.  Basic-State 

We  assume  that  each  basic  flow  quantity  is  the  sum  of  a mean  flow 
quantity  (the  Blasius  value)  and  an  unsteady  disturbance  quantity: 

U(x,y)  * Uo(y)  + e My)exp[i(kx  - wt)]  + cc  + 0(e2),  (6) 

V(x,y)  * e Vi(y)exp[i(kx  - u>t)]  + cc  + 0(e2)t  (7) 

P(x,y)  « P0  + e Pi(y)exp[1(kx  - ut)]  + cc  + 0(e2),  (8) 

where  Uo  and  P0  are  mean  flow  quantities,  c is  the  amplitude  of  the  T-S 
wave  and  cc  stands  for  the  complex  conjugate  of  the  preceding  terms.  The 
dimensionless  frequency  Is  represented  by  u and  the  real  part  of  k is 
the  wavenumber  whereas  the  Imaginary  part  of  k Is  the  negative  of  the 

growth  rate,  von  Kerczek*  treated  the  case  for  which  Eqs.  (6)  - (8)  are 
Independent  of  x. 


Substituting  Eqs.  (6)  - (8)  into  Eqs.  (1)  - (3),  subtracting  the 
mean  flow  quantities  and  linearizing  the  resulting  equations  In  the 
unsteady  disturbance  quantities,  we  obtain  the  following  equations 
describing  the  T-S  wave: 

£i(Ui,Vi;k)  s DVi  + IkUi  • 0.  (9> 

tfa(Ui,VlfPi;k,«)  = i(U0k  - u)Uj  + ViDU 0 + 1kPt  - J (D2  - k2)U,  * 0, 

(10) 

^3(U1.V1.P1;k,a>)  = 1(U0k  - u>)V,  + DP,  - \ (D2  - k2)Vi  * 0,  (11) 

. Ui  = Vi  = 0 at  y = 0,  (12) 

UlVi.Pj  -*■  0 as  y *•  «,  (13) 

where  D = 3/3y. 

B.  Stability  Analysis 

To  study  the  stability  of  the  basic  state,  we  superpose  small 
unsteady  disturbances  on  the  basic  flow  quantities  according  to 

U(x,y,t)  = U(x,y)  + u(x,y,t),  (14) 

V(x,y,t)  - V(x,y)  + v(x,y,t),  (15) 

p(x,y,t)  * P(x,y)  + p(x,y,t),  (16) 

/V  A A 

where  U,  V and  P represent  the  basic-state  given  by  Eqs.  (6)  - (8)  and 
u,  v and  p are  the  time  dependent  disturbances  which  are  assumed  to  be 
small  compared  with  the  basic-state  quantities. 

Substituting  Eqs.  (14)  - (16)  into  Eqs.  (1)  - (5),  subtracting  the 
basic-state  quantities  and  linearizing  the  resulting  equations  In  the 
unsteady  disturbance  quantities,  we  obtain 


|1  ♦ |i  . 0, 

ax  ay  * 


(17) 


,#•••  ■ 


9u  + 0^+u^+V^+v^+^-iv2u  = O, 
at  u ax  u ax  v ay  v ay  ax  r v u u* 

lv  + 0^-+u^-+V^-+v^  + ^-  -V2v  = O 
at  u ax  u ax  v ay  v ay  ay  R v v u* 


08) 

(19) 

(20) 


u = v = 0 at  y = 0, 

u,v,p  -*■  0 as  y -*•  «.  (21 ) 

To  determine  an  approximate  solution  to  Eqs.  (17)  - (21),  we  follow 
the  method  of  multiple  scales5  and  seek  a uniform  expansion  in  the  form 

u = u0(x0,Xi ,y,t)  + eui(x0,Xi,y,t)  + 0(e2),  (22) 

v = v0(x0,xi,y,t)  + evi(x0,xi,y,t)  + 0(e2),  (23) 

P = Po(x0,X! ,y,t)  + epi(xo,Xi,y,t)  + 0(e2),  (24) 

where  xi  s ex0  is  a slow  scale.  Substituting  Eqs.  (22)  - (24)  into  Eqs. 
(17)  - (21),  using  Eqs.  (6)  - (8)  and  equating  the  coefficients  of  like 
powers  of  e,  we  obtain 


Order  e° 


3ui 


(25) 


Mi(uo,v0)  = ♦ Dv0  * 0, 

M2(uo,v0,Po)  = + uo  + VoDUo  + fj£  - J (°2  - |^)v«  = °* 

(26) 

M3(uo,v0,Po)  = l^4-  + U«  v-J-  ♦ Dp0  - ^ (D2  - |^r)vo  * 0» 

u0  * Vo  5 0 at  y ■ 0, 
uo.vo.po  -*-0  as  y -*•  «, 


(27) 

(28) 
(29) 


Order  e: 

M»(ui ,Vi ) 


3U0 
3X!  • 


(30) 
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Mx(ui.vltpi)  = - Uof^J-  - [Uiexp(1kx0  - lut)  + cc]  f^J-  - 1ku0[Ui  * 
exp(lkx0  - lut)  - cc]  - [Viexp(1kx0  - 1u>t)  + cc]Du0 

- vo[OUiexp(ikx0  - Ut)  ♦ cc]  - fjfj-  ♦ £ . 

(31) 

Ns(ui.Vi.Pi)  « - u«  Uj-  - [U1exp(1kx0  - Iwt)  + cc]  - 1ku0[Vi  x 

exp(1kx0  - lut)  - cc]  - [Viexp(ikx0  - iut)  + cc]Dv0 

- v0[OViexp(1kx0  - iwt)  + cc]  + £ , (32) 

(33) 

(34) 


Ui  = Vi  = 0 at  y « 0, 
ui.vi.pi  -*■  0 as  y -*•  ». 
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III.  Solution 

A.  Basic-State  Solution 

The  solutions  to  Eqs.  (9)  - (13)  can  be  expressed  as 


Ui  = ASn(y), 

(35) 

Vi  = Aci2(y), 

(36) 

Pi  = Aci a (y ) , 

(37) 

where  A is  a constant  and  £in  (n  = 1,2,3)  are  the  eigenfunctions  of  the 
parallel  problem  given  by  the  following  equations: 

£i(tn*ti2ik)  = 0, 

(38) 

^.2(5 1 1 ,Ci2  »Ci  jjk.co)  = 0, 

(39) 

J&sUlI.ClZ.tlSik.u)  = 0, 

(40) 

Ci 1 = C12  = 0 at  y = 0, 

(41) 

Cii»Ci2»Ci3  + 0 as  y 

(42) 

B.  Zeroth-Order  Solution 


The  solution  of  the  zeroth-order  problem  given  by  Eqs.  (25)  - (29) 
is  taken  in  the  form 

uo  = B(xi )?2i(y)exp[i (Kxo  - fit)]  + cc,  (43) 

Vo  = B(xik2z(y)exp[i(Kx0  - fit)]  + cc,  (44) 

Po  * B(xi)c2j(y)exp[i(Kx0  - fit)]  + cc,  (45) 

where  K is  the  dimensionless  propagation  constant  and  fi  is  the  dimen- 
sionless frequency  for  the  secondary  mode.  The  amplitude  B is  still  an 
undetermined  function;  it  Is  determined  by  imposing  the  solvability 
condition  at  the  next  level  of  approximation.  The  quasi -parallel  Orr- 
Somnerfeld  problem  is 


l(C21,C22*.K)  » 0, 

(46) 

ot  l(C21,C22, C23jK.fi)  * 0, 

(47) 

i 
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it  j(C2l,C22,C23;K,fi)  = 0, 


(48) 

(49) 

(50) 


C21  = C22  * 0 at  y = 0, 

C21.C22.C23  "*"0  as  y -*•  °°. 

C.  First-Order  Problem 
Since  the  homogeneous  parts  of  Eqs.  ( 30 ) - ( 34 ) are  the  same  as  Eqs. 
(?5)-(29)  and  since  the  latter  equations  have  a nontrivial  solution,  the 
inhomogeneous  Eqs.  (30)-(34)  have  a solution  if,  and  only  if,  a solv- 
ability condition  is  satisfied.  The  inhomogeneous  parts  in  Eqs.  (30)- 
(34)  contain  terms  proportional  to  exp[i (kx0-  wt)],  exp[i(kxo  - wt) 

- i(Rx0  - wt],  their  complex  conjugates  and  others.  Secular  terms  will 
arise  at  this  order  when  k ~ 2K  and  w Z 2 J2;  that  is,  when  a subharmonic 
resonance  exists.  We  consider  the  case  of  perfect  time  resonance  and 
introduce  a detuning  parameter  a for  the  spatial  part  according  to 

k • 2K  eo,  0 ■ 0(1),  (51) 

w = 2 n.  (52) 

Using  Eqs.  (51)  and  (52),  we  obtain 

(kx0  - wt)  - (Rx0  - Jit)  = (Kxo  - Jit)  + (crxi  + 21K^x0),  (53) 

where  is  the  imaginary  part  and  R is  the  complex  conjugate  of  K. 

To  determine  B,  we  seek  a particular  solution  for  the  first-order 


problem  in  the  form 

ui  * Zi(y;xi)exp[i(Kx0  - Jit)]  + cc,  (54) 

vi  = z2(y;xi)exp[i(Kx0  - Jit)]  + cc,  (55) 

Pi  = z3(y;x  )exp[i(Kx0  - Jit)]  + cc.  (56) 

Substituting  Eqs.  (35)-(37),  Eqs.  (43)-(45)  and  Eqs.  ( 51 ) -(56)  Into  Eqs. 
(30)-(34)  and  equating  the  coefficients  of  exp[i(Kx0  - nt)],  we  obtain 
*£  i(zi,Z2»K)  - gi.  (57) 

£ i(zi.22,z3;K.n)  a 92.  (58) 

j(zi,z2,Zj;K,n)  * g3,  (59) 
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zi  = z2  = 0 at  y = 0 


(60) 

(61) 


zi,z2,z3  ->0  as  y + », 

where  gi,  g2  and  g3  are  defined  in  Appendix  A. 

To  obtain  the  solvability  condition,  we  multiply  Eqs.  ( 57 ) - ( 59)  by 
★ ★ ★ * 

C2 1 » c22  and  C23.  respectively,  where  the  x,  's  are  the  solutions  of  the 
adjoint  homogeneous  problem,  add  the  equations  and  integrate  the  resulting 
equation  by  parts  from  y = 0 to  y = °°.  The  adjoint  problem  corresponding 
to  the  eigenvalue  K is 


iKc22  - DC23  = 0, 

(62) 

i ( U 0 K - fi)C23  + C22DU0  - Dc23  - (D2 

- K2)ct3  = 0, 

(63) 

i(U„K  - fi)ct2+  iKc*i  - 1 (D2  - K2k*?. 

= 0, 

(64) 

C22  =5*3=0  at  y = 0, 

(65) 

* * ★ 

C2 1 2 2 1 C23  0 as  y 

(66) 

Then,  the  solvability  condition  can  be  expressed  as 

00 

f * * * 

(giC2i  + g2C22  + gaC23)dy  = 0.  (67) 

Jo 

Substituting  for  g3,  g2  and  g3  from  Appendix  A into  Eq.  (67),  we 
obtain  the  following  differential  equation  for  the  evolution  of  B: 

» ^-ABexp(ioXi  - 2^X0),  (68) 

★ 

where  fi  and  f2  are  given  in  quadratures  in  terms  of  c2n>  S2n,  k and  K 
and  they  are  defined  in  Appendix  B.  To  account  for  the  parallel  growth 
rate,  we  let 

b = Bexp(-Kjx)  (69) 

and  substitute  Eq.  (69)  into  Eq.  (68)  to  obtain 
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a^r = - Kib  + eABexp[{  ieo  - 2Ki)x0].  (70) 

So  far,  the  flow  Is  assumed  to  be  parallel.  To  account  for  the 
nonparallel  effects  only,  we  follow  Nayfeh,  Sarlc  and  Mook6  and  Sarlc 
and  Nayfeh7  and  obtain 

3x7“  (Tf  ‘ Ki)b.  (71) 

where  f3  is  defined  In  Appendix  B.  Hence,  to  account  for  the  combined 
effects  of  the  subharmonic  resonance  and  nonparallelism,  we  have 

- (ff  ‘ Ki)b  + eABexp(ieo  - 2K.)x]  (72) 

Note  that  when  the  nonparallel  effects  are  neglected  (l.e.,  fj  = 0) 

Eq.  (72)  reduces  to  Eq.  (70).  Moreover,  when  the  effects  of  subharmonic 
resonance  are  neglected,  (i.e.,  f2  = 0)  Eq.  (72)  reduces  to  Eq.  (71). 

Thus,  Eq.  (72)  describes  the  nonparallel  spatial  growth  of  the  amplitude 
of  the  subharmonic  mode  for  different  initial  amplitudes  of  the  fundamental 


1 


Two  of  these  roots  have  positive  real  parts  and  make  the  solution  grow 
exponentially  as  y -*•  « and  must  be  discarded  according  to  the  boundary 
conditions.  This  leaves  two  linear  Independent  solutions  that  decay 
exponentially  with  y.  To  use  SUPORT,  we  express  the  boundary  conditions 
at  y * y0  in  the  form 

Qz  5 0 at  y = y0»  (78) 

where  Q Is  a 2x4  matrix  with  constant  coefficients  and  given  in  Appendix 
C. 

The  eigenvalues  are  not  known  a priori  and  must  be  determined  along 

with  the  eigenfunctions.  For  given  values  of  u>  and  R,  we  guess  a value 

for  k,  generate  the  matrix  G at  y = y0  and  integrate  the  system  of 

equations  from  y = y0  to  y = 0.  If  the  guessed  value  of  k does  not 

satisfy  the  boundary  conditions  at  y = 0,  k is  incremented  by  using  a 

Newton-Raphson  scheme  and  the  procedure  is  repeated  until  the  boundary 

conditions  are  satisfied.  Integration  is  done  by  using  a technique 

8 

developed  by  Scott  and  Watts  . This  technique  orthonormal izes  the 
solution  of  the  set  of  equations  whenever  a loss  of  independence  is 
detected. 

B.  Solution  of  the  adjoint  problem 

The  solution  procedure  Is  exactly  the  same  as  for  the  solution  of 
the  basic-state  problem.  The  coefficients  of  the  z matrix  are 

zi  = t*2(y).  z2  = D;t2(y),  z3  = c!j(y),  z„  * ;2i(y),  (79) 

and  the  adjoint  problem  has  the  same  eigenvalues  as  the  zeroth-order 
problem. 

C.  Solution  of  the  solvability  condition 

The  calculations  are  repeated  at  different  x locations  to  evaluate 
fi.  f 2 » fj,  k and  K for  a given  frequency  along  the  x-axIs.  A fourth- 
order  fixed  step  size  Runge-Kutta  integration  scheme  Is  used  to  solve 
Eq.  (72)  to  find  the  amplitude  of  the  subharmonic  mode  for  different 
initial  amplitudes  of  the  T-S  mode. 
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V.  Numerical  Results  and  Concluding  Remarks 

Computations  were  performed  for  three  different  dimensionless 

frequencies  of  the  subharmonic  wave,  F = 44  x 10”®,  52  x 10”®  and 

60  x 10”®,  by  using  different  values  of  e,  where  F = ft/R  = fi*v/l£  with 

* 

ft  being  the  dimensionless  boundary-layer  frequency  and  ft  being  the 
dimensional  frequency  of  the  subharmonic  wave.  Here,  e along  with  the 
amplitude  function  A,  which  is  not  a constant  for  nonparallel  flows, 
determines  the  amplitude  of  the  fundamental  wave.  The  amplitude  function 
A can  be  calculated  by  using  the  usual  nonparallel  theory®’7  and  the 
calculated  values  of  A for  different  frequencies  are  shown  in  Fig.  1. 

If  A is  normalized  at  the  initial  location,  e then  represents  the 
initial  amplitude  of  the  fundamental  wave  as  a fraction  of  the  mean 
flow.  The  calculations  are  started  at  the  first  neutral  stability  point 
of  the  fundamental  wave  and  continued  further  downstream  well  into  the 
unstable  region  of  the  subharmonic  mode.  The  amplitude  of  the  subharmoni 
mode,  b(x),  is  also  normalized  at  the  initial  starting  point. 

Figures  2,  3 and  4 illustrate  the  variations  of  b(x)  as  functions 
of  R and  e for  three  different  values  of  F.  The  solid  curves  are  for 
the  case  of  no  interaction  between  the  fundamental  and  the  subharmonic 
waves.  The  point  where  the  solid  curve  has  a minimum  is  the  first 
neutral  stability  point  of  the  subharmonic  mode. 

The  results  show  that  the  presence  of  the  fundamental  mode  causes 
the  amplitude  of  the  subharmonic  mode  to  oscillate  before  reaching  its 
unstable  region.  As  the  subharmonic  wave  approaches  its  unstable  region, 
it  starts  to  amplify  at  a faster  rate.  Finally,  it  breaks  away  from  the 
oscillatory  growth  region  and  amplifies  considerably.  The  amplification 
of  the  subharmonic  mode  Increases  with  increasing  values  of  e. 
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The  amplification  rate  of  the  subharmonic  wave  can  be  much  faster 
than  that  of  the  fundamental  wave.  However,  as  the  amplitude  of  the 
subharmonic  builds  up,  one  can  no  longer  neglect  its  effect  on  the 
fundamental  wave  and  a nonlinear  analysis  is  necessary  to  account  for 
the  interaction. 
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APPENDIX  A 

dB 

91  = " 377  C21  (Al) 

92  = ' Uo  377  C21  " 377  523  + IT  377  521  " [i(k  " R)tu 

+ C12DC21  + Dcu?22]A6exp(ioxi  - 2^X0)  (A2) 

o - 11  dB  r + 21K  dB  r r „ - A ..  - 

93  “ " Uo  377  ^22  “R“377  ?22  ■ L-RC11C22  + ikc12£21 

+ C12DC22  + Dci2C22]A6exp(iaxi  - 2K.x0)  (A3) 


APPENDIX  B 


["  C2i£21  “ (U0S2I  + C 2 3 ) ?2  2 


" U0C22C23  + ( C2 1^2 2 


+ ?22^3)]dy  (Bl) 

^2=1  [i(k  - R)?llC21  + C12DC21  + D^ 1 1 £2 2 1 C2  2 
•'o 


+ [-  IRC11C22  + ikc12c21  + ?i2Dc22  + DCi2C22]^23  dy  (B2) 


f3  = £2  I C Czi  + C22(Uo  txf*  + Ixf  C21  + VoD?21  + 

*« 

+ (Uo  + VoDc22  + C22DVo)C23  - ^ (f^  C»  + 34f-  cJ>) 

■ 1^3x7  ^21^22  + C2aC23)]dy  (B3) 

where  x2  = e2x0  and  e2  = R"^  expressing  the  slight  nonparallelism  of  the 
flow. 
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APPENDIX  C 


9 1 1 = o,  9i2  = 1»  9 i 3 = 0,  g14  = 0 (Cl) 

921  = i(U0k  - w)R  + k2,  g22  = 0,  g23  - R , g24  ■ ikR  (C2) 
931  = - ik,  932  = g33  =934=0  (C3) 

94i  = 0,  g42  = — 1 k/R,  g43  = - [i(Ujk  - u)  + k2/R],  g44  = 0 


(C4) 

Q is  a 2x4  matrix  consisting  of  the  last  two  rows  of  the  matrix  B~^. 
The  matrix  B has  the  elements: 


bn 

= bj2  = bi 3 

= bn  = 1 

(C5) 

621 

= - k,  b22 

= k,  b23  s k,  b24  = - k 

(C6) 

bsi 

= i,  b 3 2 = 

ik/k,  b33  = - i,  b34  = - ik/k 

(C7) 

b4i 

- (w/k  - 1), 

, b42  = 0,  b43  - (oj/k  - 1),  b44  = 0 

(C8) 

where 

— 

k = 

[k2  + i(k  - 

^R]55 

(C9) 

II 


¥ 
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Figure  Captions 


Figure  1.  Amplitude  of  the  fundamental  mode  at  different  frequencies. 

Figure  2.  Amplitude  of  the  subharmonic  mode  for  different  Initial 
amplitudes  of  the  fundamental  mode  at  F = 44  x 10"®. 

Figure  3.  Amplitude  of  the  subharmonic  mode  for  different  initial 
amplitudes  of  the  fundamental  mode  at  F = 52  x 10"®. 

Figure  4.  Amplitude  of  the  subharmonic  mode  for  different  initial 
amplitudes  of  the  fundamental  mode  at  F = 60  x 10"®. 
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